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The basic geometry of a stellar interferometer with two telescopes consists of a baseline vector and 
a direction to a star. Two derived vectors are the delay vector, and the projected baseline vector 
in the plane of the wavefronts of the stellar light. The manuscript deals with the trigonometry 
of projecting the baseline further outwards onto the celestial sphere. The position angle of the 
projected baseline is defined, measured in a plane tangential to the celestial sphere, tangent point 
at the position of the star. This angle represents two orthogonal directions on the sky, differential 
star positions which are aligned with or orthogonal to the gradient of the delay recorded in the u — v 
plane. The North Celestial Pole is chosen as the reference direction of the projected baseline angle, 
adapted to the common definition of the "parallactic" angle. 
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I. SCOPE 



II. SITE GEOMETRY 



The paper describes a standard to define the plane that 
contains the baseline of a stellar interferometer and the 
star, and its intersection with the Celestial Sphere. This 
intersection is a great circle, and its section between the 
star and the point where the prolongated baseline meets 
the Celestial Sphere defines an oriented projected base- 
line. The incentive to use this projection onto the Ce- 
lestial Sphere is that this ties the baseline to sky coordi- 
nates, which defines directions in equatorial coordinates 
which are independent of the observatory's location on 
Earth. This is (i) the introduction of a polar coordi- 
nate system's polar angle in the u — v coordinate system 
associated with the visibility tables of the Optical Inter- 
ferometry Exchange Format (OIFITS) [H, [li|, detailing 
the orthogonal directions which are well and poorly re- 
solved by the interferometer, and (ii) a definition of ce- 
lestial longitudes in a coordinate system where the pivot 
point has been relocated from the North Celestial Pole 
(NCP) to the star. 

In overview, Section [TTJ introduces the standard set 
of variables in the familiar polar coordinate systems. 
Section IIIII defines and computes the baseline position 
angle — which can be done purely in equatorial or hori- 
zontal coordinates or in a hybrid way if the parallactic 
angle is used to bridge these. Some comments on relating 
distances away from the optical axis (star) in the tangent 
plane to the optical path delay conclude the manuscript 
in Section [TVl 
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A. Telescope Positions 

1. Spherical Approximation 

In a geocentric coordinate system, the Cartesian coor- 
dinates of the two telescopes of a stellar interferometer 
are related to the geographic longitudes Ai, latitudes fa 
and effective Earth radius p (sum of the Earth radius and 
altitude above sea level), 



cos Xi cos ( 
T; = p \ sin Xi cos c; 



= 1,2. 



(1) 



We add a "<?" to the geocentric coordinates to set them 
apart from other Cartesian coordinates that will be used 
further down. A great circle of radius p « 6380 km, cen- 
tered at the Earth center, connects Ti and T2; a vector 

,Ti x T 2 , where 



perpendicular to this circle is J = - 3 sinZ . 
a baseline aperture angle Z has been defined under which 
the baseline vector 



b = T 2 - Ti 

is seen from the Earth center: 



(2) 



Ti.T 2 = |T 1 ||T 2 |cosZ; (3) 
cos Z = cos <pi cos 4>2 cos( A 2 — Ai ) + sin <fii sin </> 2 .(4) 

J is the axis of rotation when Ti is moved toward T 2 . 
The nautical direction from Ti to T 2 is computed as 
follows: define a tangent plane to the Earth at TV Two 
orthonormal vectors that span this plane are 



Ni = 



- cos Ai sin<^i 

- sinAi sin</>i 

cos</>i 




dTx/dfa |Ni| = l. 



(5) 
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to the North and 



example 



Ei 




0Ti/0Ai; [Ei| - 1. (6) 



to the East. The unit vector from Tx to T 2 along the 
great circle is J x -Tx. The compass rose angle r is 

computed by its decomposition J x -Ti = cosrNi + 

sinrEx within the tangent plane. It becomes zero if T 2 
is North of Tx, and becomes n/2 if T 2 is East of Tx: Fig. 

m 




star 



FIG. 1: In the horizontal coordinate system, the baseline di- 
rection is characterized by the azimuth angle At. In the case 
shown, < A < ty < A b and D < 0. 

(Note that the tables in Q use a different definition.) 
Straight forward algebra establishes r from 

cos0i sin 02 — sin 0i cos 02 cos(Ax — A 2 ) . . 
COST = : — - ; (7) 



smr = 



sin Z 
cos 02 sin(A2 — Ax) 
sin Z 



(8) 



2. Caveats 

The angle r of the previous section does not transform 
exactly into r ± tt if the roles of the two telescopes are 
swapped, because the great circle, which has been used 
to define the direction, is not a loxodrome. This asym- 
metry within the definition indicates that a more generic 
framework to represent the geometry is useful. 

If we consider (i) geodetic rather than geocentric rep- 
resentations of geographic latitudes [l4[ or (ii) long base- 
line interferometers build into rugged landscapes, the 5 
parameters above (p, A;, 4n) are too constrained to han- 
dle the 6 degrees of freedom of two earth-fixed telescope 
positions. Nevertheless, a telescope array "platform" is 
helpful to define a zenith, a horizon, and associated co- 
ordinates like the zenith distance or the star azimuth. 
This leads to the OIFITS concept, an array center C, for 



C x 
C z 



plus local telescope coordinates 

Ti = C+t u i = 1,2. 

B. Sky coordinates 



(9) 



(10) 



Geodetic longitude A, geodetic latitude and altitude 
H above the geoid are defined with the array center 0, 

HE IE Eg! , 



C z 



(N + H) cos cos A 
(N + F)cos0sinA 
[N(l - e 2 ) +H] sin? 



(11) 



where e is the eccentricity of the Earth ellipsoid, and 
N = p e / yl — e 2 sin 2 the distance from the array center 
to the Earth axis measured along the local vertical. 

In a plane tangential to the geoid at C we define a 
star azimuth A, a zenith angle z, and a star elevation 
a = 7r/2 — z, 



- cos A sin z 
sin A sin z 
cos z 



- cos A cos a 
sin A cos a 
sin a 



(12) 



We label this coordinate system t as "topocentric," with 
the first component North, the second component West 
and the third component up. The value of A used in this 
script picks one of the (countably many) conventions: 
South means A = and West means A = +tt/2. This 
conventional definition of the horizontal means the star 
coordinates can be transformed to the equatorial param- 
eters of hour angle h = I — a, declination <5, and right 
ascension a [H (5.45)] d (2.13)], 

cos a sin A — cos<5sin/i; (13) 

cos a cos A — — sin S cos + cos 5 cos h sin 0; (14) 

sin a = sin 5 sin + cos 8 cos h cos 0; (15) 

cos 5 cos h — sin a cos + cos a cos A sin 0; (16) 

sin 6 — sin a sin — cos a cos A cos 0. (17) 

These convert (1121) into 



cos sin 5 — sin cos 5 cos h 

cos 6 sin h 
sin sin 6 + cos cos 5 cos h 

A star at hour angle h = is on the Meridian, 

sin(<5 — 0) \ 

, (h = 0), 

cos(<5 - 0) J t 



(18) 



(19) 
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north of the zenith at A = it ii S — <fi > 0, south at A = 
if S — 4> < 0. A reference point on the Celestial Sphere 
is the North Celestial Pole (NCP), given by insertion of 
5 = tt/2 into (US), 



(20) 




s+ = 

The angular distance between s and s + is 
s • s + = sin 5. 

C. Baseline: Generic Position 



(21) 



The baseline vector coordinates of the geocentric 
OIFITS system 



B x 

By 



Ti 



(22) 



(defined to stretch from Ti to T 2 — the opposite sign is 
also in use [lTj]) could be converted with 



t2x 



2y 



I =U t g{4>, A)" By 



*** 1 1 V tlz > 1 




to the topocentric system via the rotation matrix 

- sin <j> cos A — sin (j> sin A cos < 
U tg {^X)=\ sinA -cosA ]. (24) 

cos (f) cos A cos cf> sin A sin < 



Just like the star direction (|12|) . the unit vector b along 
the baseline direction defines a baseline azimuth Ab and 
a baseline elevation ab 



- cos Ab cos ab 
sin Ab cos ab 
sin a;, 



(25) 



The standard definition of the delay vector D and pro- 
jected baseline vector P is 



D 



D 



P _L D. 



(26) 



P — and later its circle segment projected on the Celestial 
Sphere inherits its direction from b such that heads and 
tails of the vectors are associated with the same portion 
of the wavefront: Fig. [2] We use the term "projected 
baseline" both ways: for the straight vector (of length P, 
units of meter) that connects the tails of D and b, or the 
line segment (length 6, units of radian) on the Celestial 
Sphere. 

The dot product of ([T2]) by $?5§ is 




D 



•b = 6cos0, (0<6»<tt) 



(27) 



FIG. 2: Sign conventions of the co-planar vectors b, D and 



where the angular distance 6 between the baseline and 
star directions is introduced as 

cos 9 = cos ab cos a cos(^4b — A) + sin ab sin a. (28) 

The star circles around the NCP in 24 hours, which 
changes the distance to the baseline in a centric periodic 
way fApp. |B|) . 



III. PROJECTED BASELINE ANGLES 
A. Definition 

We define position angles of points on the Celestial 
Sphere relative to a star's position as the bearing angle 
by which the NCP must be rotated around the star direc- 
tion s (the axis) until the NCP and the particular point 
of the object are aligned to the same direction (along a 
celestial circle) from the star. The sign convention is left- 
handed placing the head of s at the center of the Celestial 
Sphere. This is equivalent to drawing a line between the 
star and the NCP, looking at it from inside the sphere, 
and defining position angles of points as the angles in 
polar coordinates in the mathematical sign convention, 
using this line as the abscissa and the star as the center. 
Another, fully equivalent definition uses a right-handed 
turn of the object around the star until it is North of 
the star. And finally, as an aid to memory, it is also 
the nautical course of a ship sailing the outer hull of the 
Celestial Sphere, which is currently poised at the star's 
coordinates and is navigated toward the point. 

This defines values modulo 2tt; whether these are fi- 
nally represented as numbers in the interval [0, 2tt) or in 
the interval (— 7r, +tt] — as the SLAJ3EAR and the SLA_PA 
routines of the SLA library [23[ do — is largely a matter 
of taste. 

An overview on the relevant geometry is given in Fig- 
ure[3]which looks at the celestial sphere from the outside. 
In Figure [3l the baseline has been infinitely extended 
straight outwards in both directions which defines tele- 
scope coordinates Ti and T 2 also where the baseline pen- 
etrates the Celestial Sphere. Standing at the mid-point 
of the baseline, telescope 1 then is at azimuth t, telescope 
2 at azimuth Ab = r + ir, see Fig. [T] Let the object be 
at azimuth A and elevation a. 

The projection of the baseline occurs in a plane in- 
cluding the object and the baseline and thus defining the 
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FIG. 3: Celestial sphere, seen from the outside. The north 
direction through the object is given by the great circle pass- 
ing through the celestial poles and the object. The angle 8 is 
the vector P = b D projected on the sphere. 



cos A sin a 
— sin A sin a 
cos a 



(31) 



(The equivalent exercise with a different choice of axes 
follows in Section lillDI ) These are orthonormal, 



s-e A 



s-e a = e A -e a = 0; \e A \ 



1; 



sxe a = e^. 
(32) 

Decomposition of the direction from s to the NCP defines 
three projection coefficients co,i, 2 , 



s+ = c s + cie A + c 2 e Q . 



(33) 



Building the square on both sides yields the familiar for- 
mula for the sum of squares of the projected cosines, 



S, = l. 



(34) 



Two dot products of (J33J) using ([20]) and ((30J) d3TJ) solve 
for two coefficients, 

ci = s + ■ e A = cos</>sinv4, (35) 
c 2 = s + • e a = cos (/> cos j4 sin a + sin ^> cos a. (36) 



great circle labeled in Figure [3) If 9 is the angle on the 
sky between the object and the point on the horizon with 
azimuth At, then the length P of the projected baseline 
is given by 



P = b sin0, 



(29) 



where 8 is calculated from the relation (|28|) . With this 
auxiliary quantity, there are various ways of obtaining 
the position angle pb of the projected baseline on the sky, 
some of which are detailed in Sections IIII Bl MIDI The 
common theme is that the baseline orientation {Ab, a&) 
in the local horizontal polar coordinate system is trans- 
ferred to a rotated polar coordinate system with the star 
defining the new polar axis; the position angle is the dif- 
ference of azimuths between T 2 and the NCP in this 
rotated polar coordinate system. 



B. In Horizontal Coordinates 

We use a vector algebraic method to span a plane tan- 
gential to the celestial sphere; tangent point is the po- 
sition s of the star. Within this orthographic zenithal 
projection Q we consider the two directions from the 
origin toward the NCP, s + , on one hand and toward T 2 , 
b, on the other. The tangential plane is fixed by any two 
unit vectors perpendicular to the vector s of (fl"8)) . Rather 
arbitrarily they are chosen along ds/dA and ds/da, ex- 
plicitly 



e i 




(30) 




FIG. 4: In the tangent plane spanned by the unit vectors 
(|3Up and (|31l) . the components of b and s+ are ci, c^, C2 and 
c' 2 defined by (|33[) and (|40[) . The view is from the outside 
onto the plane, so e A , the axis vector into the direction of 
increasing A, points to the left. 

An angle ip between e A and s + is defined via Fig. [4] 

ci 



COS if 

sirup 



C-2 



(37) 
(38) 



5 



From (|34[) . and since (|2ip equals Co, 
\J c{ + c 2 = cos S. 



The computational strategy is to build 



(39) 



The equivalent splitting of b into components within and 
perpendicular to the tangent plane is: 

b = c s + c[eA + c' 2 e a . (40) 

Building the square on both sides yields 

Co 2 + cf + 4 2 = l. (41) 



Dot products of (|40)l solve for the expansion coefficients, 
with (ED and (I30l)-(I3H): 



c[ = b • e A 

= — cos Ab cos a;, sin A + sin Ab cos a& cos A 

= cos a b sin(^l fc - A), (42) 

c' 2 = b • e a 

= — cos a& sina cos(y4t, — A) + sinab cosa. (43) 

This defines an angle ip' in the polar coordinates of the 
tangent plane, 



COS ip 



sinyj 



Since c' equals b • s = cos 9 in 



c'2 + c', 2 = sin ( 



(44) 
(45) 

(46) 



The baseline position angle is the difference between the 
two angles, as to redefine the reference direction from ba 
to the direction of the NCP (see Fig.g]): 

Pb = (f' -<p (mod 2ir). (47) 

Sine and cosine of this imply 1, 4.3.16, 4.3.17] 



smpi, = simp cos ip — cos if simp 
c\d 2 -d x c 2 ^ 
cos 8 sin 9 ' 
cos pb = cos ip' cos tp + sin ip' sin ip 
_ cic[ + c 2 c' 2 
cos 6 sin 9 ' 



(48) 
(49) 



where 

c\c 2 — C\C 2 — cos </>(sin A sin cos a — cos ab sin a sin ^4^) 
— sin cos a cos a& sin(^4b — A), (50) 

and 

c\c x + c 2 c2 = cos0(cos A cos a cos — cos a;, cos Ab) 

+c' 2 sin(j)cosa. (51) 



tanpb = 



ClC 2 - C^C2 

cic[ + c 2 c' 2 ' 



(52) 



The (positive) values of cos 8 and sin0 in the denomi- 
nators of (|48|) and (|49|) need not to be calculated. The 
branch ambiguity of the arctan is typically handled by 
use of the atan2() functionality in the libraries of higher 
programming languages. 



C. Via Parallactic Angle 

If the position angle of the zenith p is known by any 
other sources — see App. E] — a quicker approach to the 
calculation of pb employs an auxiliary angle tp [ll| , 



Pb = p + 7r + tp (mod27r). 



(53) 



The calculation of ip is delegated to the calculation of its 
sines and cosines 

cosip = — cos pb cosp — sin pb sinp; (54) 
simp = — sin pb cos p + cos pb sin p. (55) 

In the right hand sides we insert (|48] l -(|49 |) and ([AT])- 
L2p. and after some standard manipulations 



sin^ 
cos ip 



cos ab sin(ylb — A) 
sin 9 

sin a cos 9 — sin ab 

sin 9 cos a 
cos ab sin a cos(Ab — A) — sin ab cos a 



(56) 



(57) 



These are c[ and -c' 2 of (g5]) and (TJ3]) divided by sin 9. 
In Fig. 21 ^ can therefore be identified as the angle be- 
tween — e a and b, and this is consistent with (|53[) which 
decomposes the rotation into the angle p, a rotation by 
the angle 7r (which would join the end of p with the start 
of ip in Fig. Q|, and finally a rotation by ip. 

In summary, the disadvantage of this approach is that 
p must be known by other means, and the advantage is 
that computation of ip via the arctan of ([56j) and (|57|) 
only needs c[ and c' 2l but not c\ or c 2 . 



D. In Equatorial Coordinates 



Section UlI Bl provides ph. given star coordinates (^4, a). 
Subsequently we derive the same value in terms of S and 
h. The straight forward approach is the substitution of 
CE3])-([I7]) in ((501) and §51$ . A less exhaustive alternative 
works as follows [9j : the axes and e a in the tangential 
plane are replaced by two different orthonormal direc- 
tions, better adapted to 5 and h, namely the directions 
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along along ds/dh and ds/dS: 



Building the dot product of (|63|) with e$ yields 



e/, 



sin (j> sin h 

cos h 
- cos 6 sin h 



(58) 



cos </> cos S + sin </> sin 8 cos /i 
f,j [ — sin 5 sin h | . (59) 

sin cos 8 — cos </> sin 8 cos ft. 

The further calculation follows the scheme of Section 
IIIIBI The axes are orthonormal, 

s-e h = s-es = e h -es = 0; \e h \ = \eg\ = 1; sxe5 = ed. 

(60) 

The decomposition of s+ in this system defines three ex- 
pansion coefficients do,i,2 ; 

s + = d s + d^s + d 2 e h . (61) 

Multiplying this equation in turn with and yields 

d 2 = 0; di = s + ■ e$ = cos 8. (62) 

The direction b to T2 is projected into the same tangent 
plane, defining three expansion coefficients d' 12 , 

b = d' s + d[es + d! 2 e.h- (63) 
Building the dot product of this equation with using 




FIG. 5: Fig. H] after switching from the (eA,e a ) axes to 
(eh, eg). Axis components of s+ and b are indicated. 



2B yields 

d' 2 = — cos Ah cos ah sin </> sin h + sin Ah cos ah cos h 

— sinab cos 0sin h. (64) 



d[ = — cos Ah cos ah cos 4> cos 5 

— cos Ah cos ab sin <fi sin <5 cos /i 

— sin Ah cos a;, sin 8 sin /i + sin ah sin </> cos 8 

— sina& cos ^sinJcos h. (65) 

Some simplification in this formula is obtained by trading 
the hour angle h for the angle 9, which might be more 
readily available since 9 is measured via the delay: In the 
second term we use (Ti"4l to substitute 



. , , cos 4> sin 8 + cos A cos a 

sin <p cos h — > . (66) 

cosd 

In the third term we use (Tl3|) to substitute 

. , sin A cos a , . 

sin h -> — , (67) 

coso 

and in the last term we use (fT5)) to substitute 

sin a — sin 5 sin 

cos0cosra-^ . (68) 

coso 

Further standard trigonometric identities [l], 
4.3.10,4.3.17] and (HHJ) yield 

. sin ah sin d> — cos a;, cos Ah cos — sin 8 cos 

rf i = c ■ ( 69 ) 

coso 

In the numerator we recognize a baseline declination Sh, 

s + • b = sin Sh = sin ah sin (j) — cos ah cos Ah cos 0. (70) 

is the angle that rotates the projected vector (d 2 , di) 
within the tangent plane into the direction (d' 2 , d[ ) (Fig. 



ED, 



cospb 
sinpf, 



d' 



suit 



(71) 
(72) 



Again, the sin^ does not need actually to be calculated, 
but only 



tanpt 



-4 
d[ 



(73) 



and again, selection of the correct branch of the arctan is 
easy with atan2() functions if the negative sign is kept 
attached to d' 2 . 

Note also that the use of (|69|) is optional: d' 2 and d[ 
are completely defined in terms of the baseline direction 
(Ab,ab), the geographic latitude 4> an d the star coordi- 
nates (h, 8) via l|6"4"|) and ([63)1 . With these one can imme- 
diately proceed to (|73"1) . 

OIFITS defines no associated angle in the plane 

perpendicular to the direction of the phase center. 
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The segment of the baseline in the tangent plane pri- 
marily defines an orientation in the u — v plane. However, 
modal decomposition of the amplitudes in products of 
radial and angular basis functions (akin to Zernike poly- 
nomials) means that the Fourier transform to the x — y 
plane preserves the angular basis functions (see the cal- 
culation for the 3D case, Spherical Harmonics, in [IH). 
In this sense, the angle p b is also "applicable" in the x — y 
plane. 



IV. DIFFERENTIALS IN THE FIELD-OF-VIEW 
A. Sign Convention 

The provisions of the previous section, namely 

• the OIFITS sign convention of the baseline vector, 
Eq. @ and Fig. [2J 

• the conventional formula J77|) 

fix the sign of the delay D as follows: if the wavefront 
hits T2 prior to Ti, the values of D and cos are positive; 
if it hits Ti prior to T2, both values are negative. 

The transitional case D — (6 — 7r/2) occurs if the 
star passes through the "baseline meridian" plane per- 
pendicular to the baseline; the two lines to the North- 
west and Southeast in Fig. [T] show (projections of) these 
directions. The baseline in most optical stellar inter- 
ferometers is approximately horizontal, a b w 0. From 
(|27p we see that this case is approximately equivalent to 
cos(^4fc — A) w 0. Since cos a > 0, the sign of D coin- 
cides with the sign of cos(A(, — A). This is probably the 
fastest way to obtain the sign of the Optical Path Dif- 
ference from FITS @ header keywords; the advantage of 
this recipe is that the formula is independent of which 
azimuth convention is actually in use. 

The baseline length b follows immediately from the Eu- 
clidean distance of the two telescopes entries for column 
STAXYZ of the OIFITS table 0I_ARRAY [IB]. 



B. External Path Delay 

The total differential of |28|) relates a direction 
(AA, Aa) away from the star to a change in the angular 
distance between the star and the baseline, 

— sinOAO = Aa [— cosa b sinacos(v4b — A) + sinab cos a] 
+AA cos a cos a b sin(A b - A) . (74) 

ASS]) and ([57]) turn this into 

sin Ad = sin 6>(cos ipAa — sin t/j cos a A A) , (75) 

consistent with the decomposition in Fig. [5] The two 
directions of constant delay are characterized by AO = 0. 



Equating the right hand side of the previous equation 
with zero, this means 



Ac 



cos aAA 



t&nip (A6 = AD = 0). (76) 



For the direction of maximum change in AD, which 
runs perpendicular to (|T6"|) , 



Aa 



1 



cosaAyl tani/> 
which leads to 



, ( Pb : max AD) (77) 



AD 



-bsinOAO 

, . cos aAA . 
bsmO — : — ■ — = — osm( 

smtp 



} Aa 

cos ip 



(maxAL>()78) 



Example: a field of view of AO = ±1" at a baseline 
of b = 100 m at a "random" position = 45° scans 
AD = ±340 /im. The details of the optics determine how 
far this contributes to the instrumental visibility (loss). 

Fig. [6] sketches this geometry: the change AD when re- 
pointing from one star by (AA, Aa) to another star along 
some segment of the Celestial Sphere can be split into a 
change associated with the radial direction toward/away 
from T2 along a great circle, and a zero change moving 
along a circle of radius cos centered on the baseline 
azimuth. 



zenith 




FIG. 6: Example of three concentric circles on the Celestial 
Sphere, centered at the baseline, which unite star directions of 
D = const. Two star positions are connected by the short di- 
agonal dash. A quarter of each of the two projected baselines, 
which intersect at (At,ab) near the horizon, is also shown. 

We may repeat the calculation in the equatorial system 
and split AO into AS and Ah, 

AD = -bsinOAO (79) 

= b( ASe s ■ b + Ah cos 5e h ■ b) (80) 

= b{ASd[+ Ah cos Sd' 2 ) (81) 

= 6sin6'(A(5cospb — A/icosi5sinpb). (82) 
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This translates @, App. C] to our variables. When 
AO = 0, we have the directions of zero change in AD, 
perpendicular to the projected baseline, 



AS 



tanp b , (AO = AD = 0). 



cos SAh 

The direction of maximum change is 
AS 1 



cos SAh tanp;, 
and along this gradient 

AD 



(pb : max AD), 



(83) 



(84) 



, . „cosSAh , . „ AS , . 
6sin6>— = 6sin6» , (max AD). (85) 



sin Pb 



COS Pb 



V. SUMMARY 

To standardize the nomenclature, we propose to choose 
the North Celestial Pole as the "reference" direction (di- 
rection of position angle zero) and a handedness to define 
the sign of the position angles (mathematically positive 
if looking at the Celestial Sphere from the inside) . 

The set of position angles discussed here includes 

• the direction of the zenith, the "parallactic" angle, 

• the direction toward the point where the baseline 
pinpoints the Celestial Sphere, the "projected base- 
line angle," 

• position angles of "secondary" stars in differential 
astrometry. 

Computation of the angle proceeds via (|52|) if the ob- 
ject coordinates are given in the local altitude-azimuth 
system, or via (|73p if they are given in the equatorial 
system. 

The mathematics involved is applicable to observato- 
ries at the Northern and the Southern hemisphere: posi- 
tions (A, a) or ((5, h) are mapped onto a 2D zenithal co- 
ordinate system. The position angles play the role of the 
longitude (the North Celestial Pole the role of Greenwich 
or Aries). A distance between the object and the star (in 
the range from to 7r along great circles through the Ce- 
lestial Pole in an ARC projection) might take the role of 
the polar distance — although an orthographic SIN pro- 
jection had been used for the calculations in this script. 



APPENDIX A: PARALLACTIC ANGLE 



cos a. This also turns (128D into cos 9 = sin a. (l48l) 



(gU) become 
sinp 
cosp 



i A 



cosS 

cos 6 cos A sin a + sin 6 cos a 



cos S 

The ratio of these two equations is 

cos 6 sin A 



tanp 



cos d> cos A sin a + sin d> cos a 



(Al) 
(A2) 

(A3) 



This expression can be transformed from the horizontal 
to equatorial coordinates as follows: Multiply numerator 
and denominator by cos a and use (|13[) in the numerator, 



tanp = 



cos 6 cos S sin h 



cos (j) cos A sin a cos a + sin <j> cos 2 a 
Insert (|17[) in the denominator, 

cos <t> cos S sin h 

tanp 



(A4) 



sin a(sin a sin <f> — sin 6) 
cos 4> cos 8 sin h 
sin d> — sin a sin S 



sin m cos^ a 



(A5) 



Replace sin a with (JT5J) in the denominator, eventually 
divide numerator and denominator through cos</>cos<5: 



tanp 



cos 6 cos 8 sin h 



sin — (sin 5 sin 4> + cos S cos h cos 4>) sin S 

cos (j> cos S sin h 
cos 2 S sin <f> — cos S sin S cos h cos <j> 

sin h 



cos S tan d> — sin S cos h 



(A6) 



These manipulations involved only multiplications with 
positive factors like cos 0, cos S and cos a; therefore one 
can use sin h and cos S tan <fi — sin 5 cos h separately as 
the two arguments of the atan2 function to resolve the 
branch ambiguity of the inverse tangent. 

This definition of the parallactic angle coincides with 
the dominant one in the literature U [TH, H3, H2l j 
others exist: Equation (1) in Q generates the 90°- 
complement of our definition, for example. 

If we insert ([57 ) -((56] ) into (JTSJ) [1, 4.3.18] and compare 
with (SI, we find 



AS 



cos aAA 



tanp 



cos 5 Ah 1-tanp 



(A7) 



cos aAA 



which means the parallactic angle mediates between di- 
rections expressed as axis ratios in the two coordinate 
systems (12| . 



We define the parallactic angle p as the position angle 
of the zenith. Since the formulas in Section HTT1 dealt with 
the general case of coordinates on the Celestial Sphere, 
one may just set a,b = n/2 in (|25p. which induces 0^ = 0, 



APPENDIX B: DIURNAL MOTION 

The diurnal motion of the external path difference is 
described by separating the terms oc sin/i and cx cosh. 
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We rewrite (1271) in equatorial coordinates by multiplying 
([15]) with (J2SJ, and collect terms with the aid of ([70]) : 

D = s • b = 6[sin5sin<5b + cos 5 cos 5b cos (h — hb)]. (Bl) 

A time-independent offset b sin <5 sin 5b and an amplitude 
& cos <5 cos (*>(,, defined by products of the polar and equa- 
torial components of star and baseline, respectively, con- 
stitute the regular motion of the delay. The phase lag hb 
is determined by the arctan of 

cos 5b sin hb — sin Ab cos ab; (B2) 
cos 5b cos hb — cos Ab cos ab sin0 + sina& cos 0,(B3) 

and plays the role of an interferometric hour angle. 

APPENDIX C: CLOSURE RELATIONS 

An array Ti , T2 and T3 of telescopes forms a baseline 
triangle 

bi2 + b 23 + b 31 = 0. (Cf) 

If all point to a common direction s, some phase closure 
sum rules result: 

D 12 +D 2 3 + D 3 i =0; (C2) 
D u + D 23 + D 31 = 0; (C3) 
612 cos #12 + 623 cos 023 + 631 cos 31 = 0. (C4) 

The triangle of projected baselines in the tangent plane 
is 

P12 + P23 + P31 = 0. (C5) 

These vectors can be represented as numbers in a com- 
plex plane by introduction of three moduli P and three 
orientation angles pb'- 

Pi 2 e' Pb12 + P 23 e' Pb2:i + P 3 ie ipi > 31 

= bn sin^xae'f" 2 + 6 23 sin 9 23 e vp ^ 3 + b 31 sin9 31 e lpb31 

= 0. (C6) 

These equations remain correct if transformed by com- 
plex conjugation or multiplied by a complex phase fa.ctor: 



therefore these closure relations of the projected baseline 
angles are correct for (i) both senses of defining their 
orientation, and (ii) any reference direction of the zero 
angle. 



APPENDIX D: NOTATIONS 



a right ascension 

a, Aa star elevation and its difference 

A, AA star azimuth angle, and its difference 

Ab baseline azimuth 

b, b baseline vector and length 
5 declination 

D signed external path difference 

eh, eg, e a, e a unit vector in tangent plane 

e eccentricity of geoid 

4> geographic or geodetic latitude 

h hour angle 

H altitude above geoid 

I local sidereal time 

A geographic longitude 

N distance to the Earth axis measured 

along local normal 
p position angle of the zenith, "parallactic 

angle" 

Pb position angle of the projected baseline 

P projected baseline length 

p Earth radius 

p e Earth equatorial radius 

s unit vector into star direction 

s + unit vector into NCP direction 

9 angular separation between star and 

baseline vector 
tp angle between projected baseline circle 

and star meridian 
z zenith angle 

Z baseline angle seen from the Earth cen- 

ter 
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